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LONG-TERM  GOAL 

The  long-tenn  goal  is  to  develop  a  new  powerful  capability,  which  is  called  SNOW  (simulation  of 
nonlinear  ocean  wave-field),  for  predicting  the  evolution  of  large-scale  nonlinear  ocean  wavefields 
using  direct  phase-resolved  simulations.  Unlike  the  phase-averaged  approaches,  SNOW  models  the 
key  physical  mechanisms  such  as  nonlinear  wave-wave,  wave-current,  wave-wind  and  wave-bottom 
interactions  and  wave-breaking  dissipation  in  a  direct  physics-based  context. 

OBJECTIVES 

The  specific  scientific  and  technical  objectives  are  to: 

1 .  Develop  effective  physics-based  modeling  of  wind  forcing  and  wave-breaking  dissipation  for 
robust  phase-resolved  computation  of  nonlinear  wavefield  evolution. 

2.  Extend  SNOW  capabilities  to  handle  high  sea  states  for  investigating  the  effect  of  very  steep  local 
waves  upon  evolution  of  wave  statistics  of  nonlinear  wavefields. 

3.  Extend  SNOW  to  general  finite  water  depth  by  including  effects  of  bottom  dissipation,  fluid 
stratification,  and  variable  current  and  bottom  topography. 

4.  Speed  up  the  computational  algorithm  underlying  SNOW  simulations  for  large  spatial-temporal 
scale  wavefields 
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5.  Obtain  direct  validation  and  quantitative  cross-calibration  of  SNOW  simulations  with  phase- 
averaged  wave  model  predictions  and  field/laboratory  measurements. 

APPROACH 

SNOW  employs  direct  physics-based  phase-resolved  simulations  for  predicting  the  evolution  of  large- 
scale  nonlinear  ocean  wavefields.  SNOW  is  fundamentally  different  from  the  existing  phase-averaged 
models  in  that,  under  SNOW,  key  physical  mechanisms  such  as  wave-current,  wave-wind  and  wave- 
bottom  interactions  and  wave-breaking  dissipation  are  modeled,  evaluated  and  calibrated  in  a  direct 
physics-based  context.  In  SNOW,  detailed  phase-resolved  information  about  the  wavefield  is 
obtained,  from  which  the  statistical  wave  properties  can  also  be  derived. 

SNOW  is  based  on  an  extremely  efficient  high-order  spectral  (HOS)  approach  for  direct  computation 
of  nonlinear  ocean  wavefield  evolution.  HOS  is  a  pseudo-spectral-based  method  that  employs 
Zakharov  equation  and  mode-coupling  idea  and  accounts  for  nonlinear  wave-wave,  wave-current,  and 
wave-bottom  interactions  to  an  arbitrary  high  order  (M)  in  wave/bottom  steepness.  This  method 
obtains  exponential  convergence  and  (approximately)  linear  computational  effort  with  respect  to  M  and 
the  number  of  spectral  wave/bottom  modes  ( N ).  SNOW  is  an  ideal  tool  for  phase-resolved  prediction 
of  realistic  ocean  wavefield  evolution. 

By  incorporating  point  and/or  whole  field  wave  measurements  into  the  simulations,  SNOW  provides  a 
capability  of  reconstructing  and  forecasting  nonlinear  evolution  of  phase-resolved  ocean  wavefields. 
The  objective  of  wave  reconstruction  is  to  obtain  detailed  specifications  (including  phase)  of  a 
nonlinear  wavefield,  which  matches  given  (directly  or  remotely  sensed)  sensed  wave  data  or  specified 
wave  spectrum.  Nonlinear  wave  reconstruction  is  achieved  based  on  the  use  of  optimizations  with 
multiple-level  (theoretical  and  computational)  modeling  of  nonlinear  wave  dynamics.  Using  the 
reconstructed  wavefield  as  initial  conditions,  SNOW  simulation  would  provide  a  deterministic 
forecasting  of  the  phase-resolved  wavefield  evolution.  The  validity  of  the  nonlinear  reconstruction  and 
forecasting  methodology  has  been  systematically  verified  against  laboratory  measurements  and 
synthetic  wave  data  for  both  long-  and  short-crested  irregular  wave-fields  (Wu  2004;  Yue  2008;  Wu  et 
al  2009a,  b). 

SNOW  computations  can  now  be  routinely  performed  for  nonlinear  ocean  wavefields  in  an  domain  of 
0(10  ~  )  km"  with  an  evolution  time  of  0(1)  hours.  Such  large-scale  SNOW  simulations  are  normally 
performed  on  advanced  high-performance  computing  platfonns  using  up  to  O(103)  processors  (Wu, 

Liu  &  Yue  2005;  Xiao,  Liu  &  Yue  2009)  under  our  DoD  challenge  project:  “Large-Scale  Deterministic 
Predictions  of  Nonlinear  Ocean  Wavefields”. 

WORK  COMPLETED 

The  effort  of  the  research  is  mainly  focused  on  (i)  development  of  physics-based  wind  input  modeling 
in  the  SNOW  simulation,  (ii)  extension  of  SNOW  to  account  for  effects  of  high  sea  states,  variable 
current,  density  stratification,  and  bottom  topography,  (iii)  speedup  of  SNOW  computations  on 
advanced  high-performance  computing  platfonns,  (iv)  validation  of  SNOW  simulations  by  direct 
comparisons  against  laboratory  experiments  and  field  measurement,  and  (v)  application  of  SNOW 
simulations  to  understand  generation  mechanisms  and  characteristics  of  rogue  waves.  Specifically, 
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•  Modeling  of  wind  forcing  input  We  continue  the  development  and  validation  of  wind  input 
modeling  for  SNOW  simulations  of  nonlinear  ocean  wavefield  evolution.  A  phenomenological 
wind  input  model  based  on  an  external  pressure  forcing  on  ocean  surface  is  developed  and 
incorporated  into  SNOW  simulations.  The  pressure  distribution  is  parameterized  according  to  the 
sheltering  theory  and  comparisons  to  the  laboratory  and  field  measurements  of  the  spectral  growth 
rates.  The  effectiveness  and  limitations  of  this  model  are  investigated  by  quantitative  comparisons 
with  available  measurements. 

•  Development  of  an  efficient  algorithm  for  steep  waves:  A  highly  efficient  computational 
algorithm,  so-called  pre-corrected  FFT  method  (PFFT),  is  developed  and  applied  for  the 
simulation  of  fully-nonlinear  steep  wave  dynamics.  This  approach  is  based  on  the  boundary- 
element  method  with  the  use  of  the  fast  Fourier  transfonn  technique  to  accelerate  the  evaluation  of 
influence  coefficients  (Yan,  Liu  &  Yue  2006).  With  PFFT,  the  requisite  computational  effort  in 
solving  the  nonlinear  boundary  value  problem  is  reduced  from  O (N  ~  )  to  O (N  In  N),  similar  to 
that  in  HOS.  Significantly,  PFFT  allows  the  simulation  of  fully-nonlinear  wave  dynamics.  This  is 
a  useful  complementary  to  SNOW.  Integration  of  this  algorithm  into  SNOW  would  extend  the 
capability  of  SNOW  for  fully-nonlinear  simulations  of  extreme  wave  dynamics  and  high  sea  states 
while  retain  the  high  computational  efficiency  retained. 

•  Investigation  of  stratified  fluid  and  bottom  topography  effects  upon  wavefield  evolution:  We 
extend  and  apply  SNOW  simulations  to  littoral  zones  including  stratified  fluid  and  bottom 
topography  effects.  To  consider  the  density  stratification  effect,  SNOW  is  extended  to  multi-layer 
fluids.  In  particular,  the  resonant  interactions  among  surface  waves,  interfacial  waves,  and  bottom 
ripples  are  extensively  investigated.  The  study  provides  an  understanding  of  alternate 
mechanisms  for  the  generation  of  internal  waves  in  the  ocean,  and  establishes  a  framework  for 
large-scale  phase-resolved  computations  of  internal  wavefield  evolution  and  interaction  with 
surface  waves  (Alam,  Liu  &  Yue  2009a,  b).  In  addition,  we  also  study  the  high-order  resonant 
interactions  of  three-dimensional  surface  waves  with  bottom  ripples,  which  help  understand  the 
complexity  of  wavefields  in  littoral  zones,  in  particular,  the  generation  of  infragravity  waves 
(Alam,  Liu  &  Yue  2009c). 

•  Speedup  and  applications  of  SNOW  simulations:  We  constantly  improve  the  computational 
speed,  scalability  and  robustness  of  the  SNOW  code  on  HPC  platforms  for  the  simulation  of  large- 
scale  ocean  wavefield  evolutions.  We  continue  to  apply  large-scale  SNOW  computations  to 
investigate  the  characteristics  of  statistic  quantities  of  nonlinear  ocean  waves  and  to  understand 
the  generation  mechanisms  and  statistical  features  of  rogue  waves  in  deep  ocean. 

RESULTS 

The  main  results  include:  (i)  improvement  and  perfonnance  of  SNOW  modeling  of  key  physical 
processes  such  as  wind  input,  wave  breaking  dissipation,  and  nonlinear  wave  interactions  with  density 
stratification  and  bottom  topography;  (ii)  understanding  of  nonlinear  (triad  and  quartet)  resonant 
interactions  among  surface  waves,  internal  waves,  and  bottom  topography  for  the  prediction  of 
nonlinear  surface/intemal  wavefield  evolution  in  littoral  zones  (Alam,  Liu  &  Yue  2009a,  b,  c);  (iii) 
further  validations  of  SNOW  simulations  of  nonlinear  wave  spectrum  evolution  by  comparisons  to 
phase-averaged  model  predictions  (Wu  et  al  2009);  and  (iv)  understanding  of  nonlinear  wave  grouping 
dynamics  in  the  development  and  occurrence  statistics  of  rogue  waves  (Xiao  et  al  2009).  We  below 
describe  sample  results  on  the  perfonnance  of  wind  input  modeling  as  well  as  the  effects  of  resonant 
wave-bottom  interactions  upon  ocean  wavefield  evolution  in  more  details. 
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SNOW  simulation  of  nonlinear  ocean  wavefield  evolution  with  wind  forcing:  With  the  surface 
pressure  based  model,  SNOW  simulation  can  effectively  include  the  wind  input  energy  to  the  ocean 
wavefield  and  properly  capture  the  key  characteristics  of  nonlinear  wind-sea  evolution.  In  particular, 
SNOW  simulation  provides  a  confirmation  of  wind-sea  equilibrium  hypothesis  and  Toba’s  3/2  law. 

Figure  1  shows  a  sample  SNOW  simulation  result  on  the  time  variation  of  the  total  energy  and  peak 
period  of  the  wavefield  during  the  evolution  of  a  long-crested  wavefield  under  the  influence  of  a  wind. 
The  initial  wavefield  is  given  by  a  JOHNSWAP  spectrum  with  the  effective  steepness  kpArms  =  0.107. 
At  time  t/Tp  =  500,  a  wind  forcing  with  wave  age  Cp/u*=10  and  forcing  level  akp/gpw  =  0.0004  is 
turned  on,  where  a  is  the  ratio  between  the  surface  pressure  and  the  wave  slope  in  wind  modeling. 

In  the  initial  evolution  before  the  wind  forcing  is  turned  on,  only  weak  and  isolated  wave  breaking 
occurs.  Energy  of  the  wavefield  decays  slightly  and  the  downshift  of  the  peak  wavenumber  obtains. 
After  the  wind  forcing  is  added,  in  the  initial  growth  phase  (t/Tp  =  500  to  -750),  the  wavefield  absorbs 
all  of  the  wind  energy  input  with  a  smooth  growth  in  the  wave  amplitude.  In  this  phase,  no  significant 
increase  in  the  frequency  of  wave  breaking  is  observed  and  the  growth  of  the  wave  amplitude  follows 
closely  the  prediction  by  the  linear  generation  theory  (as  indicated  by  the  green  curve  in  figure  la). 
After  the  initial  growth  phase  (t/T  r  >  -750),  wave  breaking  becomes  intermittent  with  gradually 
increasing  strength  and  frequency.  As  figures  la  and  lb  show,  the  wavefield  reaches  wind-sea 
equilibrium  with  both  wave  energy  and  peak  wavelength  (or  period)  continuing  to  grow  as  most  of 
wind  energy  input  is  balanced  by  the  wave  breaking  dissipation.  The  SNOW  simulation  result  shows 
that  in  this  continuous  growth  phase,  the  significant  wave  height  and  peak  period  of  the  wavefield 
follow  the  relation  given  by  Toba’s  3/2  law. 

Oblique  sub-  and  super-harmonic  Bragg  resonance  of  surface  waves  by  bottom  ripples:  Oblique 
quartet  Bragg  resonance  interaction  of  surface  waves  with  bottom  topography  is  considered  and  a 
select  of  possible  scenarios  are  investigated  via  theoretical  analysis  and  direct  simulation  underscoring 
their  importance  in  real  ocean  application.  At  the  third  order  in  perturbation  expansion,  two  free 
surface  waves  can  be  in  resonance  via  two  bottom  components  (Class  II  Bragg  resonance),  or,  three 
surface  waves  via  one  bottom  component  (Class  III  Bragg  resonance).  Contrary  to  Class  I  and  II, 
where  participating  waves  in  a  triad/quartet  resonance  have  the  same  frequency,  in  Class  III  they  are 
sub/super-harmonic  of  incident  wave  frequencies.  Therefore  Class  III  can  offer  a  variety  of 
possibilities  for  the  generation  of  new  frequency  resonant  waves,  hence,  affecting  the  evolution  of 
ocean  wave  spectrum. 

Specifically  it  is  shown  that  energy  can  transfer  from  short  waves  to  ordinary  ocean  waves  via  typical 
nearshore  sandbars  and  in  the  presence  of  reflected  waves  the  similar  sandbars  can  resonate  shore 
nonnal  infragravity  waves.  If  the  incident  spectrum  is  oblique,  long  longshore  waves  are  formed  that 
upon  further  class  III  interaction  can  result  in  shore-normal  semi-standing  waves.  Finally  a  case  of 
multiple  resonances  is  studied  where  simultaneous  energy  exchange  results  in  a  number  of  new  longer 
waves;  not  expected  if  class  III  Bragg  resonance  is  not  taken  into  account. 

Figure  2  presents  a  sample  result  of  direct  simulations  of  high-order  Bragg  resonant  interactions  of 
three-dimensional  waves  with  bottom  undulations.  In  this  example,  successive  third-order  Bragg 
resonances  involving  multiple  wave  and  bottom  ripple  components  lead  to  the  development  of 
irregularity  of  surface  wave  patterns  in  the  near  shore  areas. 
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Three  incident  waves  (with  wavenumers  k  1,2,3 h  =  1.0,  0.3,  0.5,  water  depth  h  =  8  m)  travel  obliquely 
(incident  angle  On  =  0.57  rad)  over  shore-parallel  bottom  ripples  with  three  components  (kbi.b2.b3h  = 
1.06,  0.72,  0.34).  Bragg  resonant  wave-bottom  interactions  generate  three  waves  (wavenumers  kri , 
kr2,  and  kr3)  traveling  in  different  directions.  Strong  reflected  waves,  with  (wavenumber,  direction) 
given  by  (krih  =  0.60,  0ri  =  2.5  rad)  and  (kr2h  =  0.40,  0r2  =  2.4  rad),  are  generated  by  third-order 
Bragg  interactions  of  three  wave/ripple  components  (ki ,  k2,  kbi)  and  (ki ,  k3,  kb2),  respectively.  After 
kri ,  kr2  waves  are  developed,  Bragg  resonant  interaction  of  kri ,  kr2  and  kb3  waves/ripples  generates  a 
transmitted  wave  with  (kr3h  =  0.20,  0r3  =  0.58),  which  is  (almost)  parallel  to  incident  k  1,2,3  waves. 

The  chain  of  multiple  resonances  continues  and  results  in  a  complex  wave  field.  Figure  2a  shows  the 
time  variation  of  the  amplitudes  of  the  resonance  generated  waves  ari ,  ar2 ,  and  ar3 .  The  amplitudes  of 
these  waves  can  reach  ~40%  of  the  incident  wave  amplitude  after  an  evolution  of  t/T  i~  50.  Figure  2b 
shows  the  complex  combined  wavefield  at  t/T  1  =50. 

IMPACT/APPLICATIONS 

This  work  is  the  first  step  toward  the  development  of  a  new  generation  of  wave  prediction  tool  using 
direct  phase-resolved  simulations.  It  augments  the  phase-averaged  models  in  the  near  term  and  may 
serve  as  an  alternative  for  wave-field  prediction  in  the  foreseeable  future. 
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Figure  1.  Phase-resolved  SNOW  simulation  of  nonlinear  evolution  of  a  long-crested  wavefield  in 
the  presence  of  a  wind  forcing:  (a)  variation  of  the  energy  of  the  wavefield  (normalized  by  the  initial 
value)  and  (b)  variation  of  normalized  peak  wavenumber  of  the  wavefield  as  a  function  of  time 
during  the  evolution.  The  wind  forcing  is  turned  on  at  t/Tp  =  500.  The  green  line  in  (a)  represents 
the  exponential  growth  of  wavefield  predicted  by  the  linear  wind-wave  generation  theory. 
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Figure  2:  Multiple  oblique  high-order  Bragg  resonances  of  waves  with  bottom  undulations  obtained 
from  direct  simulation.  Three  incident  waves  (wavenumers  k i,  k2,  and  k  {)  travel  obliquely  over 
shore-parallel  bottom  ripples,  and  Bragg  resonant  wave-bottom  interactions  generate  three  waves 
(wavenumers  krl,  kr2,  and  kri)  traveling  in  different  directions,  (a)  Growth  of  the  (normalized) 
amplitudes  of  the  Bragg-resonance  generated  waves;  and  (b)  a  snapshot  of  the  complex  wave 
pattern  resulted  from  the  combination  of  the  incident  and 
resonance-generated  waves  at  time  t/T  /  =  50. 
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